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ABSTRACT 

Using H I absorption alone, we attempt to separate H i absorption lines in quasar 
spectra into two categories; HDLs (Higher Density Lines) and LDLs (Lower Density 
Lines), and we discuss the difference in their physical properties. We deblend and fit 
all H I lines with Voigt profiles, and make an unbiased sample of H i lines covering a 
wide column density range (12 < log N^j < 19 cm~^). To reduce the influence of line 
blending, we simultaneously fit several Lyman series lines. As a result of a two-point 
correlation analysis, we found that higher column density H i lines are clustering at 
Av < 200 km s~^, while lower ones at Av < 100 km s~^. We define HDLs as H i lines 
with 15 < log Nhi < 19 cm~^ and all H i lines within ±200 km s~^ of a line with 
log Nhi > 15 cm~^, and LDLs as others with 12 < log Nhi < 15 cm~^. We found that 
the HDLs have smaller minimum 6- values for a given column density than the LDLs. 
This difference is successfully reproduced by our Hydrodynamic simulation. The LDLs 
seem to be cool or shock-heated diffuse IGM gas, while the HDLs are likely to be cooler 
dense gas near to galaxies. 
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1. Introduction 

Quasars have been used as background sources, allowing the study of objects that lie between 
us and them, most of which make Lya absorption lines. To date, line-profile fitting analysis for 
Lya absorption lines has been carried out to study physical parameters such as the velocity width 
and the column density (Pettini et al. 1990; Carswell et al. 1991; Rauch ct al. 1992; Hu et al. 1995 
(H95 hereafter); Lu et al. 1996a (L96 hereafter); Kirkman & Tytler 1997 (KT97 hereafter)). 

Samples of H i absorption lines usually contain (i) H I lines originating in the intergalactic 
diffuse gas, and (ii) H I lines produced near or in intervening galaxies. The former could be weak 
H I lines, widely distributed as the Lya forest (hereafter we call these "LDLs" or "Lower Density 
Lines"), whereas the latter are likely to be strong H i lines clustered as metal absorption lines 
(hereafter we call these "HDLs" or "Higher Density Lines"). HDLs are thought to be caused by 
discrete clouds. On the other hand, LDLs are probably arising from photoionized IGMs in the 
continuous density fields that are broadened by Hubble flow (e.g., Rauch 1998; Kim et al. 2002). 
Here, the reader can think of density as column density Nhi, or volume density, since numerical 
simulations show that the two are correlated (Zhang et al. 1998). The HDLs arc lines formed in 
or near to higher density parts of the IGM. As the log Nhj increases, especially above about 17, 
these regions are more likely to be in or near the outer parts of galaxies, or clumps that will become 
galaxies. In contrast, the LDLs arise in relatively lower density regions of the IGM. 

There are studies which suggest that some H I absorption lines have a close relationship with 
galaxies near the line of sight, not only on small scales at law-z (Grogin & Geller 1998; Penton, 
Stocke, Sz Shull 2002) but also on large scales at high- 2; (Adelberger et al. 2003). McDonald, 
Miralda-Escude, & Cen (2002) presented theoretical predictions of the correlation between the Lya 
forest transmitted flux and the mass of absorbers within ~ 5h^^ Mpc (comoving) of the line of sight. 
Other spectroscopic and imaging observations have been carried out in order to study the statistical 
relationship between the Lya absorbers and galaxies at law-z. The most important result is that 
the Lya equivalent widths are anti-correlated with the projected distance of the nearest galaxies 
within 500h~^ kpc of the galaxies (Lanzctta ct al. 1995; Tripp, Lu, & Savage 1998; Dave et al. 
1999; Chen et al. 2001). Thus, some H I absorbers seems to be closely related with galaxies. 

Thus, the analysis of HDLs and LDLs is indispensable for the detailed investigation of physical 

properties of H i gas. However, one of the most serious problems in such an analysis is that the 
Lya absorption lines are so heavily blended with each other that it is difficult to separate and fit 
them individually. This problem is often seen in strong H I lines with log Nhi > 15 cm^^. 

In some past studies, the column densities of strong H I lines have been evaluated using the 
Lyman limit optical depth for LLSs, or the total rest-frame equivalent width for DLA systems 
without fitting them by Voigt profiles (e.g. Lanzetta et al. 1991; Petitjean et al. 1993). These 
methods, however, evaluate only the total column densities and Doppler parameters of the heavily 
blended H I lines. 
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In this study, we attempt to fit all H I lines with Voigt profiles, construct an unbiased sample 
of H I lines over a wide column density range (12 < logA^^^-j < 19 cm~^), and investigate their 
physical parameters, such as column density, Doppler parameter, and clustering properties. In order 
to separate the heavily blended (strong) H i lines and fit them individually with Voigt profiles, we 
use not only the Lya line but also higher Lyman series lines, such as Ly/? and Ly7, to improve the 
fitting accuracy. 

This method has previously been applied to only a several LLS and DLA systems (e.g. Songaila 
et al. 1994; Tytler, Fan, & Buries 1996; Wampler, Baldwin, & Carswell 1996; Carswell et al. 1996; 
Songaila, Wampler, & Cowie, 1997; Buries & Tytler 1998a, 1998b, Buries, Kirkman, & Tytler 1999; 
Kirkman ct al. 2000; O'Meara et al. 2001; Kirkman et al. 2003) and the Lya forest in two quasar 
spectra (Kim ct al. 2002). Here we apply this fitting method to 40 quasar spectra acquired with 
KECK + HIRES (Vogt et al. 1994). 

Finally, we attempt to separate H I lines into HDLs and LDLs, and apply our statistical analysis 
to these two classes separately. We confirm that the column density distribution and the Doppler 
parameter distribution of LDLs are similar to the past results (H95; L96; KT97). We found the 
most remarkable difference between HDLs and LDLs in the plot of column density vs. Doppler 
parameter, which suggests that H l absorbers are not produced by a single phase (or a single 
population). This difference is also reproduced by our Hydrodynamic simulation. The summary of 
the other results and the detailed description of each absorption system are presented in Misawa 
et al. (2004). 

We give a brief description of the observation in § 2. In § 3, we explain how to evaluate the 
line parameters. In § 4, we investigate the physical properties of HDLs and LDLs, and compare 
the observational results with the Hydrodynamic simulation. We summarize our results in § 5. 

2. Observations 

The sample quasars were originally selected from a survey taken for measurements of the D/H 
(deuterium to hydrogen) abundance ratio. The typical ratio of D/H is so small, 3 xl0~^ (Kirkman 
et al. 2003 and references therein), that we can only detect D I lines corresponding to H I lines with 
large column densities, log Nhi > 16.5 cm~^. Given this fact, we observed quasars in which either 
DLA systems or LLSs were detected. Here we examine KECK+HIRES spectra of 40 quasars. We 
used a 1".14 slit which provided a resolution of 8.0 km s~^. The spectra were extracted by the 
automated program, MAKEE package, written by Tom Barlow. 
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3. Preparation of uniform sample 

Instead of detecting all the H i lines in the spectra of the 40 quasars, this study includes only 
the H I lines with log Nhi > 15 cm~^, and all other H I lines within it 1000 km of such H I 
lines. The reason is described below. 



3.1. Detection of H I systems 

We at first searched for H I lines with log Nhi > 15 cm~^, in the following way: (1) we found H I 
lines already discovered in DLA systems or LLSs in earlier work (Sargent, Steidel, &; Boksenberg 
1989; Lanzctta 1991; Tytlcr 1982; Buries 1997), (2) we checked for H i lines corresponding to 
previously discovered metal absorption systems (Peroux et al. 2001; Storric-Lombardi et al. 1996; 
Petitjean, Ranch, & Carswell 1994; Lu et al. 1993; Steidel & Sargent 1992; Lanzetta et al. 1991; 
Barthel, Tytler, & Thomson 1990; Steidel 1990; Sargent et al. 1980; Sargent, Steidel, & Boksenberg 
1988), and (3) we finally attempted to detect H I lines ourselves that satisfied both of the following 
conditions — (i) the H I line has a column density of log Nhj > 15 cm~^, and (ii) at least one 
higher Lyman series, e.g. Ly/3 and Ly7 lines, is present along with Lya in the observed wavelength 
range, which improves the reliability of line fitting. We do not consider metal lines. 

Metal absorption lines (Sargent et al. 1980; Young, Sargent, & Boksenberg 1982; Sargent, 
Steidel, & Boksenberg 1988; Petitjean & Bergeron 1994; Churchill & Vogt 2001; Pichon et al. 
2003) including those in DLA systems (Lu et al. 1996b) cluster over Av < 400 km s~^. It is then 
sufficient to detect all lines within ± 1000 km of the H I lines with log Nhi > 15 cm~^ for the 
purpose of studying the physical properties of weak and strong H I lines. Where more than one H I 
line with log Nhi > 15 cm~^ was detected within this velocity window, the position of the line with 
the largest column density (hereafter the "main H I component") was considered to be the center 
of the H I system. As an example, we show the velocity distribution of H I lines for the system at 
Zabs= 3.321 in the spectrum of Q0014+8118 in Figure 1 (see also Buries, Kirkman, & Tytler 1999). 
With this method, we detected 86 H I systems at 2.1 < z < 4.0 in the spectra of 31 quasars, out of 
a total of 40 examined. 

3.2. Evaluation of Line Parameter 

For the fitting of each absorption line with a Voigt profile, we fit all the accessible lines in the 
Lyman series, helping us deblend H I lines and make an unbiased sample of H I lines. 

In the process of line detection, we removed very narrow lines with Doppler parameter of 
b < 4.81 km s~^ or FWHM ~ 8.0 km s~^, which is the resolution of our spectra. We decided to 
ignore all lines with 6 < 15 km s~^ as described later. If the Voigt fitting allowed different solutions 
during the fitting trials, the result with the minimum number of lines was chosen. This occurs 
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mainly in the saturated region for which we have only small number of Lyman series information. 
If we fit the line profiles using smaller number of components than the actual number, the 6- values 
could be overestimated. On the other hand, the column density cannot be easily overestimated, 
because only a small rise of the column density significantly changes the line profiles (e.g., damping 
profile) especially at the saturated regions. 

Once the fitting model was applied to the absorption lines in some region of the spectra, we 
used minimization to find the model parameters which best fitted the observed spectrum. The 
formal errors on the fitted parameters for lines with 6 < 30 km s"^ are a{logNHi) = 0.09 cm~^, 
a{b) = 2.1 km s~^ and a{z) = 2.5 x 10~^. Before beginning our statistical analysis, we prepared 
an unbiased sample of H I lines, by eliminating inappropriate H I systems that have disadvantages 
described below. 

(1) poor fitting due to gaps in echelle formatted spectra: Echelle-for matted spectra some- 

times have data gaps between echelle orders. This occurs frequently in the redder part of 
the spectrum (e.g., A > 5000 A in our data). Spectral gaps appear every ~ 100 A, and the 
widths of the wider gaps are no less than 20 A depending on the wavelength. If the Lyman 
series line damaged by the spectral gaps is Lya, the fitting accuracy for the system is very 
low. Therefore we removed H I systems whose Lya lines are affected by spectral gaps wider 
than 100km s~^. 

(2) poor fitting due to strong DLA wings: DLA (and sub-DLA) systems with large H I col- 

umn densities of log Nfji > 19 cm~^ are unsuitable for our study, not only because their 
absorption profiles are damped so strongly that almost all weak components are blanketed 
by the damping wings, but also because the spectra around DLA lines are not normalized 
correctly due to the strong absorption features. 

(3) close proximity in redshift to the background quasars: The Lya forest disappears in 

the regions redward of the Lya emission lines of the quasars. Therefore we remove H I 
systems at a distance of < 1000 km from the Lya emission lines, because the asymmetrical 
distribution of their H I lines complicates the clustering analysis described in § 4.1. 

(4) overlapping with other H I systems: H I system pairs, whose velocity windows of ±1000 

km s~^ are overlap, are also excluded, as their distribution of H I lines could affect each other 
and the line clustering analysis is contaminated. 

Finally, we have an unbiased sample of 973 H I lines in 61 H I systems. However the sample is 
not homogeneous because the S/N ratio varies. The S/N ratios of the spectra are at least S/N ~ 
11 per 2.1 km s~^ pixel and the mean value is S/N ~ 47 for Lya lines. 
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4. Physical properties of H I lines 

We used the sample prepared in § 3, to study the H i hue parameters. Our analysis is similar 
to that of previous studies (e.g. H95; L96; KT97), but with three key differences: (i) earlier studies 
used all H l lines detected in the quasar spectra, whereas we use only H I lines within it 1000 km s""^ 
of the main components with log Nhi > 15 cm~^, (ii) our sample contains a number of strong H I 
lines {log Nhi > 15 cm~^) in addition to weak lines (log Nhi < 15 cm'^), and (iii) our sample 
covers a wide redshift range, 2.0 < z < 4.0, compared with those of the past studies, Az ~ 0.5. 



4.1. Clustering Properties 

Our sample contains not only LDLs but also HDLs. Dave et al. (1999) noted in their hydro- 
dynamic simulations that galaxies tend to lie near the dense regions that are responsible for strong 
H I lines. Therefore, we attempted to classify the lines using their clustering along the line of sight 
by constructing the two-point correlation function (Sargent et al. 1980), 

where N{v) is the number of observed pairs at a velocity separation of v, and Nexp{v) is the number 
of line pairs expected if they are randomly placed along the line of sight. Nf.j.p{v) is calculated by 
Monte Carlo simulations. Webb (1987) found significant clustering, ^{v) = 0.32 ± 0.08 over the 
velocity range 50 < f < 150 km s~^ at 1.9 < z < 2.8. H95 confirmed the same trend with ^{v) 
= 0.17 lb 0.045, over the same velocity range at {z) = 2.8. With a large sample of H i lines with 
log Nhi > 13.8 cm~^ at 1.7 < 2; < 4.0, Cristiani et al. (1997) also found that the correlation 
function at v ~ 100 km s~^ increases with log Nhi- On the other hand. Ranch et al. (1992), L96, 
and KT97 did not find any clustering for similar velocity ranges at 2.7 < z < 3.4, [z) = 3.7, and 
(z) = 2.7, respectively. 

In Figure 2, we show the number of H I lines of our sample at a given distance from the main 
components. The H I line with relatively large column densities tend to cluster around the main 
components, and have a symmetrical distribution, suggesting that they are related to each other. 
In contrast, the number of H I lines with smaller column densities decrease near the center of H I 
systems. 

Therefore we calculated ^(w) for some subsamplcs, varying the column density ranges. At first, 
a number of artificial lines matching the number of observed lines were randomly inserted into the 
61 H I system windows. Since the observed H I systems always have main components at v = 
km s"^, we also placed an artificial line at the center of each system in order to produce unbiased 
simulated data. This process was repeated 250 times to obtain an average value for Nexp{v). 

As a result of the analysis using all H I lines, we found that ^{v) shows a slight number excess 
at Av < 200 km s~^. The values of ^{v) for H I lines with various column density ranges in bins 
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with a spacing of 50 km s-^ at 50 < u < 200 km s"^ are summarized in Table 1. Column (1) 
is the range of column density, columns (2), (3), and (4) are the values of ^(f) and la Poisson 
errors for velocity separations of 50 - 100, 100 - 150, and 150 - 200 km s^^, respectively, and 
column (5) is the velocity width where the lower la deviation of ^(u) first goes below ^{v) = over 
V > 50 km s~^. Column (6) is the number of H I lines in each subsample. We see non-zero ^{v) at 
velocity separations of v ^ 100 km s"^ for weak H I lines, and v ~ 200 km for strong H I lines. 
The correlation degree at v = 50 - 100 km s~^ was found to have a maximum for H I lines with 
15 < log Nhi < 19 cm~^. This is one reason why we chose to use log Nhi = 15 cm~^ as part of 
our definition of HDLs. 

Within ~0.5/i~^ comoving Mpc of Lyman Break Galaxies (LBGs), the IGM contains less 
neutral hydrogen gas compared with the cosmologically averaged value at z ~ 3 (Adelberger et 
al. 2003). Some cosmological hydrodynamic simulations (e.g., Croft et al. 2002; Kollmeier et al. 
2003) showed that this effect could be produced not by UV flux from the LBGs but by the galactic 
winds. This decreased absorption was detected in Adelberger et al. (2003) because the lines of 
sight to the quasars and the LBGs are different. In our spectra, strong H i lines (maybe produced 
in the intervening galaxies) absorb most of the quasar flux within 50 km s^^ of their centers, which 
prevents us from detecting H I deficit in the vicinity of the strong H i lines. 



4.2. HDLs and LDLs 

We classified HDLs and LDLs with the following procedure. At first, we regarded all H i lines 
with log Nffi > 15 cm~^ as HDLs, because the line clustering is stronger as log Nhi increases as 
Cristiani et al. (1997) found, and we found above. But we emphasize that the value of log Nhi 
= 15 cm~^ is not strict. These strong H I lines may also be accompanied by weak H I lines that 
are physically associated with them. In fact, metal absorption lines sometimes have a core-halo 
structure; the strongest line is at the center of the absorption system, while the weak lines exist 
almost symmetrically on both sides of the strong one (e.g. Lu et al. 1996b; Procliaska et al. 2001; 
Misawa et al. 2003). Therefore, we defined a velocity distribution width for HDLs, vhdl- Since 
H I lines with relatively large column densities of log Nhi = 13 - 19, 14 - 19, 15 - 19, and 16 - 19 
cm~^ are correlated with each other within the velocity separation of ~ 200 km s"~^, we adopted 
Vhdl = ±200 km s~^, and regard all H I lines with 12 < log Nhi < 15 cm~^, within ± 200 km s~^ 
of H I lines with log Nhi > 15 cm~^ as HDLs. Hereafter we call the lines with log Nhi < 15, lying 
within 200 km s^^ of a log Nhi > 15 line, halo HDLs in distinction from the original HDLs. Finally 
all the 973 H i lines were separated into 306 HDLs (including halo HDLs) and 667 LDLs based on 
these criteria. 
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4.3. Column density — Doppler parameter relation 

Using the subsamples of HDLs and LDLs, we investigated the correlation between the Doppler 
parameters and column densities. Our results are summarized in Figure 3. 

In this figure, we also show lines with 5 < 6 < 15 km s~^. We find 62 LDLs and 16 HDLs with 
12 < log Nhi < 14 cm-2 and 5 < 6 < 15 km s'^, or 8% of the 973 LDLs + HDLs with 6 > 15 
km s~^. The proportion of HDLs to LDLs with 6 < 15 km s"^ (16/62 = 25.8%) is similar to that 
for lines with 6 > 15 and the same column density range (142/583 = 24.4%). Since we did not 
identify metal lines, many of them could be metals, and since our spectra have a range of S/N, 
many could be erroneous. Past studies that did consider these two effects found very few H I lines 
with these parameters (KT97; H95; L96). Hence, we ignore all lines with 6 < 15 km s~^. 

When we ignore H l lines with b < 15 km s~^, we found a positive correlation between the 
minimum Doppler parameter, brain: ^^(^ column density not only for LDLs approximately fitted by 



bmin = 4.0 X log 



Nhi (cm- 



+ 14.0 km s-^ (2) 



1012.5 

at 12.5 < log Nhi < 15.0 cm~^ as seen in the past papers (e.g. KT97), but also for HDLs fitted by 

'Nhi (cm-2)" 



bmin = 1.3 X log 



1012.5 



+ 10.5 km s"^ (3) 



at 12.5 < log Nhi < 19.0 cm~^. At log Nhi > 15 cm~^, the HDLs often have b values below the 
extrapolation of eqn. (2) for the LDLs. We, however, cannot compare the bmin values of LDLs 
with those of HDLs directly at log Nhi > 15 cm~^, because there are no LDLs at that region. 
Nonetheless, if we look at the region of log Nhi = 14 — 15 cm~^, a difference in the distribution of 
LDLs and HDLs is clear. At log Nhi = 14 - 15, only 2 % (2 of 91) LDLs are below the eqn. (2), 
while 20 % (12 of 59) HDLs are located there. We do not think this difference is just a statistical 
accident. At log Nhi < 14 cm~^, we do not see any remarkable difference between them, which 
could mean that there is little difference between the physical properties of LDLs and halo HDLs. 



4.4. Comparison with simulation 

For weak H I lines, the correlation between column density and minimum Doppler parameter 
has already been reproduced by CDM simulations. Zhang, Anninos, Sz Norman (1995) performed 
the hierarchical three-dimensional numerical simulation, and found the minimum value is increas- 
ing slightly with column density in a linear fashion. In order to mimic observational analysis more 
closely, Zhang et al. (1997) synthesized absorption spectra, and extracted, deblend, and fit the 
absorption features in generated spectra. The recovered data exhibit a bmin dependence on the 
column density well fitted by bmin = 5.5 x log[A?^(cm^^)/10^^-^] -|- 12.8 km s~^, which is consistent 
with the KT97's and our results for LDLs (eqn. (2)) shown in Figure 3. 
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We have performed a three-dimensional hydrodynamical simulation of Lya clouds in a CDM 
dominated universe, and compared it with our observational results. We used the following input 
parameters; rif,=0.04, ^7^=0. 30, Q\=0.70, Hq=70 km s^^ Mpc^-*^, (78 =0.73, and power spectrum 
slope n = 1. We chose the size of the computational box to be L= 0.7h~^ Mpc with the effective grid 
resolution of 256^. We used the Haardt &; Madau (1996) quasar spectrum with the photo-ionization 
rate, jhi = 5.6 x 10~^^ x (1 + z)^-^^ x exp[— (z — 2.3)^/1.95] ionizations per H I atom per second in 
the optical thin limit. Self-shielding is not included in the simulation. Figure 4 shows the number 
of H I lines at a given distance from the main component that are detected in our Hydrodynamic 
simulation at z = 2.9 — 3.0 (similar figure to Figure 2). The distribution of H I lines are examined 
for three column density ranges, 13.5 < log Nhj < 14.5 cm~^, 14.5 < log Nhi < 15.5 cm~^, and 
15.5 < log Nffi < 16.5 cm~^. H I lines with relatively large column density, log Nhi > 14.5, tend 
to cluster around the main components. On the other hand, weaker H I lines show the decrease in 
their number near the main components. These trends are similar to those seen in the observation 
results, although the components with 14.5 < log Nfjj < 15.5 cm~^ seem to be more concentrated 
in velocity in the simulated spectra. 

In Figure 5, we show scatter plots between log Nhi and b from the simulation. The difference of 
the distributions between HDLs and LDLs that we found in the observational result are successfully 
reproduced; strong HDLs with log Nhi > 15 cm~^ have Doppler parameters smaller than the 
extrapolation of b^nin 

for LDLs (eqn. (2)). 

One of the differences between the observation and the simulation is that there is a blank region 
at log Nhi > 15 cm~^ and 6 > 30 km only in the simulation. It may be because absorption 
lines at the region in the observation are not completely resolved even with our improved fitting 
method (especially in the case that only a few orders of Lyman series are available) , and have large 
6- values compared to the values that we might see with very high S/N ratio. The minimum 6- values, 
however, may be evaluated correctly even in such cases. On the other hand, the failure of the line 
deblending rarely happen in the simulation, because we have not included noise in the synthetic 
spectra, and because the resolution of our simulation, 256^ 0.7h^ Mpc (Ax = 2.7 kpc), is high 
enough compared to the minimum resolution (Ax = 37.5 kpc) necessary to evaluate appropriate 
6- values (Bryan et al. 1999). 

Another difference is that few narrow HDLs, whose Doppler parameter is smaller than the 
value given by eqn. (2), are detected at log Nhi < 15 cm~^ in the simulation (hereafter weak- 
narrow HDLs). Here, we would like to emphasize that the 6- values of H I lines near the strong 
H I absorbers could be over-estimated in the simulation since we did not include self-shielding and 
shadowing effects that can be effective for the HDLs. We discuss this in the next section. 
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5. Summary and Discussion 

In this paper, we, for the first time, separated H i absorption hnes into HDLs and LDLs, 
and we find differences between them using 40 quasar spectra. For weak H I lines (i.e. LDLs), 
the correlation between column density and the minimum Doppler parameter have been often 
investigated in both of the observation and the simulation. It is thought that this relation arises 
by the mechanism that higher column density clouds appear to be associated with denser regions 
in which the gas is adiabatically compressed and heated (e.g. Kim, Cristiani, & D'Odorico 2001, 
2002). Kim, Cristiani, & D'Odorico (2001, 2002) also suggested that for weak H i lines the slope 
of the log Nh I — bjnin relation became flatter as z decreased at 2; > 3, though there was no such 
significant trend at lower redshift 2; ~ 2. We, however, do not find any redshift evolution in our 
sample (Figure 6 and 7). 

The log Nhi — bmin relation for HDLs is clearly shallower than the extrapolated relations for 
weak H l lines in KT97 and Zhang et al. (1997), while the relation for LDLs is good agreed with the 
past results. This means that HDL absorbers have relatively small Doppler parameters compared 
with the LDL absorbers at the same column density. 

Dave et al. (1999) classify the H l gas clouds into three phases; (i) a cool low-density phase, (ii) 
a shock-heated intermediate-density phase, and (iii) a cold dense phase in galaxies. The first phase 
with relatively low column densities of log Nhj < 15 cm~^ make the positive correlation between 
log Nh J and bmin, while the last phase with large column densities makes an anti-correlation. 

When we think about H i clouds in the cold dense phase, the effect of radiative cooling must 
be taken into consideration. Such dense clouds with log Nhi > 17 cm~^ (i.e. LLSs and DLA 
systems) can shield themselves against the background UV flux, which preserves more H I (i.e. 
shielding effect). A similar effect may happen in the low column density clouds if they are located 
near the high column density clouds (i.e. shadowing effect). These effects arc important at least at 
z > 6 (Nakamoto, Umemura, & Susa 2001), although they have not yet evaluated quantitatively 
at lower redshift because of high complexity after the formations of first stars and quasars. Zhang 
et al. (1995) also pointed out the importance of the shielding effect, though they did not include 
the effect in their simulation, and nor did we here. If plenty of H I gas is preserved, the cooling 
of H I gas becomes very effective in regions in photoionization equilibrium (H I line cooling is the 
most effective cooling factor at T ~ 10"^ K), which could lower the temperature, and decrease the 
value of the minimum Doppler parameter. Both of these effects will tend to produce weak-narrow 
HDLs in the log Nhi — i> plane, which makes the distribution of H i lines in the simulation more 
consistent with the observation for the HDLs. 

We would like to acknowledge to M. Chiba, R. Nishi, I. Murakami, M. Nagashima, and K. 
Okoshi who gave us valuable comments on the theoretical interpretations. This work was supported 
in part by NASA grant NAG5-13113. We wish to thank the anonymous referee for the report, which 
improved the clarity of this presentation. 
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Fig. 1. — Observed and modeled velocity map for 10 Lyman series lines of the H I at Zq^^ — 3.321 
in the spectrum of Q0014+8118. The lower line in each panel is the la error. Tick marks denote 
the positions of H I absorption lines, and the large one at the center is the position of the main 
component. 

Fig. 2. — Distribution of 973 H I lines, including the main components, in 100 km s"""*^ wide bins 
from the main components. 

Fig. 3. — Doppler parameter vs. column density. Open and filled circles are LDLs and HDLs, 

respectively. LDLs and HDLs with 5 < 15 km s~^ are marked with small circles. Solid line is the 
relation between the column density log Nhi and the minimum Doppler parameter 6mm in KT97 
(12 < log Nhi < 15 cm~^), and dotted line is the same relation extrapolated up to log Nhi = 18 
cm~^. Dashed line is the same relation for HDLs in our sample. The pure thermal Doppler 
parameters corresponding to the temperatures, T = 10^, 5 x 10^, 10^, and 2 x 10^ K are shown 
with the thin dot-dashed horizontal lines. 

Fig. 4. — Same as Figure 2, but the result from the Hydrodynamic simulation, except that the 
column density ranges of three subsamples are slightly different. 

Fig. 5. — Similar to Figure 3, but from the Hydrodynamic simulation. Crosses are weak H i lines 
with log Nhi < 15 cm""^. Open square and triangle are completely deblended (6 < 40 km s~^) H I 
lines with the column densities of 15 < log Nhi < 17 cm~^ and log Nhi > 17 cm~^, respectively. In 
contrast with Figure 3, We do not use the open squares for all lines within 200 km s~^ of lines with 
log Nhi > 15 cm^^. Solid-and-dotted line and dashed line denote the relation between log Nhi 
and bjnin for LDLs (eqn. (2)) and HDLs (eqn. (3)) from the observational results. 

Fig. 6. — Column density vs. Doppler parameter relation for HDLs at 2; > 2.9 and z < 2.9. We 
chose z = 2.9 when dividing H I sample into two subsamples because it makes the numbers of H I 
systems in them nearly identical. 

Fig. 7. — Same as Figure 6, but the result for LDLs. 



Table 1. Clustering properties of H I lines 



(1) (2) (3) (4) (5) (6) 
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^correlation and the la Poisson error for the velocity separations of 50 - 
100, 100 - 150, and 150 - 200 km s'^. 

''velocity width at which the lower la deviation of ^{v) first goes below 
^{v)=0 over -y > 50 km s~^. 

^number of H I lines. 

'^we prefer to list v = 200 km s~^, though the lower la deviation of ^{v) 
goes below ^{v)=0 at 50 < < 100 km s~^; this could be due to the small 
sample. 
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Fig. 1 : Velocity map of H I system at Zabs — 

3.321 in the spectrum of Q0014+8118. 
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Fig. 2 : Distribution of H I lines in 100 km s~^ bins from the main components. 
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Fig. 3 : Column density — Doppler parameter relation for the observation. 
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Fig. 4 : Same figure as Figure 2, but the result from the simulation. 
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Fig. 6 : Column density — Doppler parameter relation for HDLs at z > 2.9 and z < 2.9. 
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Fig. 7 : Same figure as Figure 3, but the result for LDLs. 



